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Tunelowanie elektronów w złączach metal—izolator— metal 
Туннелирование электронов в системе М—И—М 


The Electron Tunneling in MIM Junction 


WSTĘP 


W większości prac dotyczących problemu tunelowania w złączach me- 
tal—izolator—metal, metal opisywano modelem Sommerfelda, natomiast 
izolator traktowano jako prostokątną barierę potencjalną. Zgodność wy- 
ników doświadczalnych z teorią bazującą na powyższym modelu uzy- 
skiwano wprowadzając parametry, takie jak „grubość efektywna” war- 
stwy izolatora oraz „masa efektywna” tunelujących elektronów. Podejś- 
cie takie nie dawało możliwości dokładnego zrozumienia zjawiska tune- 
lowania. 


5.1. Kurtin i współprac. [1] zaproponowali, aby przy interpretacji 
wyników doświadczalnych posłużyć się współczynnikiem transmisji elek- 
tronów przez warstwę izolatora, obliczonym przy pomocy metody Zene- 
ra [4]. W pracy [2] obliczono numerycznie współczynnik przenikalności 
D dla złącza, w którym izolator był opisany potencjałem Kroniga-Pen- 
neya. Prostym, ale bardziej realnym modelem potencjału jest poten- 
сја} kształtu cost-cosy-cosz. 

Celem tej pracy jest znalezienie współczynnika D dla złącza o takim 
kształcie potencjału, zarówno w metalu, jak i w izolatorze, i porówna- 
nie go ze współczynnikiem Zenera obliczonym dla tego samego modelu 
złącza. 
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MODEL 


Przyjmujemy potencjał krystaliczny postaci (ryc. 1). 


Ryc. 1. Model potencjału krystalicznego 
złącza użyty w pracy 


V(x) =V u) [8 (—x)+0 (x—z,)]--V,(x)6 (х—х,) 0 (x) (1) 
gdzie 


П 
Ум(<)=2 Ум (i сов 2 z) — Eo 


M 


П 
V,(1)=2W, (1 +сов2— z) 
І 
ам, а — stałe sieci metalu (M) i izolatora (I), 4Wy, 4%; — amplitudy 
wahań potencjału (1) w odpowiednich obszarach, Ө(х) — funkcja schod- 
kowa, znikająca dla x < 0. 


WSPÓŁCZYNNIK TUNELOWANIA 


Współczynnik tunelowania D obliczamy jako stosunek prądów: prze- 
chodzącego przez izolator ję do padającego jpaa. Operator prądu jest zde- 
finiowany następująco: 


вр h dy 
= = R ZZA э 
j=j де | im dz (2) 
gdzie w jest funkcją spełniającą jednowymiarowe równanie Schródin- 
gera oraz warunki brzegowe: a) funkcja y i jej pochodne są ciągłe w 
punktach nieciągłości potencjału, b) w obszarze III зр reprezentuje falę 
biegnącą od bariery. Równanie Schródingera w obszarach I, II, III przyj- 


muje postać równania Mathieu: 
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d?y(z 
à РЕ + [4,—, 2q, cos 221] y (2) = 0 (3) 
gdzie: ż=1,2 odpowiednio dla metalu M — obszar I, III oraz izolatora 


I — obszar II. 


2 та? 


E=E, EB=0 


4 та? 
а= PII? W, 
П 
Zu” аге Z~ Fy SE 


Rozwiązanie tego równania przedstawiamy w postaci [3] 
-+20 +еә 
(2) = Ае! w „APR еі" +к)2 + Ве > = Бачае kot (4) 


Т = —оо T=- 

W ostatnim wzorze k=0 lub 1, zależnie od parzystości odpowiedniego 
obszaru rozwiązań równania Mathieu [3]. Stałe р, Cə+x są funkcjami 
energii Е i parametrów a, W, W obszarze rozwiązań niestabilnych 
(przerwy energetyczne w krysztale) należy w rozwiązaniu (3) zastąpić 
В przez іи gdzie м jest rzeczywiste, dodatnie. Współczynniki Сз: staną 
się oczywiście liczbami zespolonymi. Po wstawieniu funkcji (4) do rów- 
nania Mathieu (3) znajdujemy związki rekurencyjne pozwalające wyzna- 
czyć parametry р, и oraz współczynniki Сук. 

W przypadku tunelowania elektronów przez pasmo dozwolone izola- 
tora otrzymujemy na współczynnik D wyrażenie: 


it 113 = 
D(E)= [coso е Ё + га sin? o| (5) 
4 р: 

_ м DFB) Cr к 

B am У С.к У, ФТ ЕВ, к 


П 
adr 


natomiast w przypadku tunelowania przez obszar zabroniony izolatora 


—2 


| (6) 


Р(Е) = icoshp + 5 Wad sinh 


24 — Annales 
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a, D (2rt+k+P)C; к Уа, ук 


0м. УС > Еа, . 
П 
ф= и — Х 
еш. 


-.W powyższych wzorach Dbaręk, Qar4k OZNACZAJĄ Współczynniki rozwinięcia 
funkcji Mathieu w obszarach energetycznie dozwolonych i niedozwolo- 
nych w izolatorze, natomiast сә: współczynniki rozwinięcia w metalu. 


METODA ZENERA OBLICZANIA D 


Stosowane obecnie metody przybliżone obliczania współczynnika D 
są metodami WKB modyfikowanymi w ten sposób, że wektor falowy 
k=V2m (E-V(x))/h jest zastępowany wektorem falowym k(E) wyzna- 
czonym z rozwiązania równania Schródingera dla izolatora. Zależność 
E(k) w okszarach dozwolonych energetycznie jest wyznaczana w stan- 
dartowy sposób (struktura pasmowa). W pasmach zabronionych do- 
puszcza się rozwiązania o wektorze falowym zespolonym, a część uro- 
joną wektora Imk(E) interpretujemy jako współczynnik tłumienia [2]. 

Współczynnik Zenera [4] tunelowania międzypasmowego jest dany 
wzorem: 


Ta 

DE)=exp| — a| Imk (Gaz (7) 
х, 

gdzie L=x— x, jest długością drogi tunelowania. W naszym przypadku 

wartość k(E) w obszarze dozwolonych energii jest rzeczywista, zatem: 


Рр(Е) = 1 (8) 


W obszarze zabronionym w izolatorze тату: k(E)=iull/a, oraz L=a. 
Korzystając z faktu, że k nie zależy od « znajdujemy ze wzoru (7) 


D(E)= exp| — 2 4 (9) 
І 


WYNIKI 


W celu zbadania dokładności metody Zenera obliczono numerycznie 
współczynnik tunelowania dla przyjętego modelu w oparciu o dokładny 
wzór (5) oraz na podstawie przybliżonego (7). Obliczenia wykonano dla 
tunelowania elektronów przez pierwsze i drugie dozwolone pasmo izo- 
latora. 
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Do obliczeń numerycznych przyjęto następujące parametry 

ам=3 А, а=6 А, Wy=0.25eV, W,=1.36eV, Eọ=2 eV 
Parametr 8 oraz współczynniki С».к występujące we wzorze -па D li- 
czone były rachunkiem samouzgodnionym z dokładnością do 1074. Ra- 
chunki przeprowadzono dla trzech różnych szerokości izolatora 24 A, 
48 A, 198 A. Zależność dokładnego współczynnika D od energii oraz jego 
zmiany z grubością warstwy przedstawia ryc. 2. W rozważanych obsza- 
rach energii (pasma dozwolone) współczynnik Zenera ma wartość D=1. 


Ryc. 2. Wykres zależności współ- 
czynnika tunelowania od energii; 
poszczególne krzywe odpowiadają 
szerokościom izolatora: a — 24А, 

b — 48 A, c — 198 A 1,9 22 М, 53 58 6,3 Е (eV) 


Widać wyraźne różnice między współczynnikiem obliczonym dokład- 
nie, a współczynnikiem Zenera. Z tego powodu należy z dużą ostroż- 
nością odnosić się do wyników uzyskanych metodami tunelowania z wy- 
korzystaniem współczynnika Zenera. 
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PE3IOME 


В настоящей работе дана энергетическая зависимость вероятности 
прохождения электронов через изолятор в системе металл—изолятор— 
металл (M—M—M). Кристаллический потенциал системы принят в BM- 
де потенциала Матье. Найденный в этой работе коэффициент проник- 
новения электронов сравнивается с коэффициентом вытекающем из 
теории Зенера. 


БОММАВУ 


The energy dependence of the exact analytically calculated transmis- 
sion coefficient for the electron tunneling in metal—insulator—metal 
(МТМ) junction is given. The transmission coefficient calculated for Ma- 
thieu crystal potential model is compared with that obtained from the 
approximation method — commonly used in the interpretation of the 
experimental data. 


